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Abstract 

The dd aX reaction is studied in a model where each pair of 
nucleons in the projectile and target deuterons undergoes pion pro- 
duction through the NN — > dvr reaction. The condition that the two 
deuterons fuse to form an a-particle then leads to peaks at small 
missing masses, the well-known ABC enhancement, but also a broad 
structure around the maximum missing mass. With a simplified in- 
put amplitude the model gives a quantitative description of both the 
a-particle momentum and angular distributions for a deuteron beam 
energy of 1250 MeV. 
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Measurements of the a-particle momentum spectra in the dd — > aX re- 
action reveal rich structure for beam energies throughout the 0.8 — 1.4 GeV 
region [Q. In particular there are very sharp peaks near the kinematic limits, 
i.e. the maximum in the a-particle cm. momentum, and these correspond 
to missing masses of Mx ~ 300 MeV. In addition, however, there is a broad 
central bump around the maximal missing mass, i.e. in the region where 
the a-particle is at rest in the cm. system. Rather similar structures were 
observed many years ago in both the pd '^HeX and np ^ dX 0, ^ 
reactions and this low missing-mass peak is commonly known as the ABC- 
effect. It is clear from the mass of the ABC that it must consist of two pions 
and the absence of an equivalent peak in the pd —>■ ^HX' case [0] shows that 
it must have isospin 1 = 0. It is therefore likely that the ABC manifests itself 
so clearly in deuteron-deuteron collisions because it leads to a pure 1 = 
channel for the pions. 

The authors of ref. originally suggested that the ABC might corre- 
spond to a large s-wave isoscalar tt-tt scattering length but phase shifts ob- 
tained from high energy pion production experiments show this is in fact very 
small p. Since the mass and width of the ABC peak vary with kinematic 
conditions, it is natural to see if it can be understood as a kinematic effect 
associated with the presence of nucleons or nuclei. In this spirit Risser and 
Shuster explored the np — ^ d X reaction in a model where two different A- 
resonances were excited through pion exchange. Though they could generate 
peaks in the right regions, the calculations had limited quantitative success, 
but it is not clear whether this was due to the simplifications in their formal- 
ism, and in particular their neglect of spin 0]. Certainly double p-wave 
pion production should lead to kinematic enhancements in the spinless case 
because the simplest matrix element squared is proportional to (fci ■ ^2)^, i.e. 
the square of the scalar product of the two pion momenta. This is maximal 
when the two momenta are parallel (the ABC peak) and antiparallel (the 
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central bump). 

ABC production in deuteron-deuteron scattering is at its largest around 
600 MeV per nucleon and this is close to the maximum of single-pion pro- 
duction in pp dTT"*" arising from A-excitation. It therefore seems natural 
to consider two-pion production as coming from two such independent reac- 
tions, involving two different pairs of nucleons from the projectile and target 
deuterons, where the a-particle is formed when the two final deuterons stick 
together. If we neglect the Fermi motion in the initial deuterons, the cm. sys- 
tems of the deuteron-deuteron and nucleon-nucleon channels coincide so that 
the produced deuterons (and pions) will have the same cm. momenta. The 
relative momentum in the final d-d system will therefore be very small when 
the pions come out with low excitation energy in the ABC peak and there will 
then be little suppression coming from the a:dd sticking factor. This clas- 
sical argument of favourable kinematics provides another reason why ABC 
production might be seen clearest in the dd aX case. 

The Feynman diagram corresponding to our model for dd aTT'^n' is 
shown in fig. |l|, where the various momenta are defined in the cm. system. In 
the non-relativistic limit of the Fermi momenta, the matrix element deduced 
from this diagram is 

'y'^ [271)2 J 

" i:^(-^')AC^^M| + (AB,«Ai5.). (1) 

spins V ^ 1 ' 

where the (e, e') are the polarization vectors of the internal deuterons and 
(ei, £2) those of the incident. The wave functions of the final a-particle and 
the initial deuterons ipaiQa) &^d (pd{q'i) depend upon the Fermi momenta 
in the rest systems of the corresponding nuclei. In the deuteron case the 
longitudinal component is boosted with respect to the cm. system to give 
Qi = {.Q\iQ.ih)i where 7 = E^/yiia and h and z denote the transversal and 
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longitudinal components. The analogous relativistic effect is small for the a- 
particle and so qa = Qi — Q2~ (^i ^ ^2)/2. Since our model is not dependent 
upon exotic Fermi momentum components, it is reasonable to neglect the 
D-state components in the deuterons and the a-particle. 

The matrix element in eq. (|l|) is essentially that of second order per- 
turbation theory, where the energy denominators corresponding to the two 
different time orderings are 

/\Ei = T {vA-{qi + q2) + {uJi-uJ2)/2} , (2) 

with va being the cm. deuteron velocity and Ui the pion total energies. This 
simple result only follows when retaining just the linear terms in the nuclear 
Fermi momenta. In such a case the two principal value integrals in eq.(|I|) 
cancel to leave only the 5-function term. 

The largest pp d7r+ amplitude in the A region corresponds to the ^D2P 
transition for which 

=C [3(p-fc)(p-e^) - fc-e^] , (3) 

where p and k are the proton and pion cm. momenta and C is a function 
of energy. This form does not depend upon the proton spin variables and 
it is due to this that we have the simplification in eq.(0) which leaves only 
an ei • €2 dependence. Keeping only this term, the value of C is fixed by 
experimental pp — > d7r+ data in the forward direction, as summarised in the 
SAID database g]. 

Some care must be taken over the effect of the a-particle binding energy 
in determining the energy in the proton-proton system. At the ABC-peak the 
pions emerge with zero relative momentum and we assume that the labora- 
tory kinetic energy in the 7r"*"d system in the inverse reaction is one half of that 
for the tt+tt^q; when the two pions have equal momentum vectors. With this 
prescription the interval of deuteron kinetic energies Td = 787 — 1412 MeV 
is mapped onto a proton range of the pp d7r+ input Tp = 410 — 720 MeV. 



Since in the linear Fermi momentum approximation only the (5-function 
survives the integration of eq.(|l]), we define the form factor 

(4) 

This integral is most easily evaluated by transforming into configuration 
space: 



W = — ^ I bdb <lzidz2 $d(6, 7^1) '^'d(-&, 7^2) ^1 (b, 



X 



Jo{\kl - k'^\b/2} cos{{kl - kl){zi - ^2)/4} cos{(a;i - u:2){zi + Z2)/Av^} . (5) 

After averaging over the spin directions of the initial deuterons, the matrix 
element squared becomes 

1 ^2 7^2 

9 E l-^r = -ff {3(p-fc0(p-fc2) + fci-fc^}' , (6) 

where A^^ is the number of deuteron pairs in the a-particle. In the ABC 
region where ki ^ k2, this yields a characteristic (3cos^6'q, + 1)^ angular 
distribution, which is just the square of that for the pp dn^ differential 
cross section in the case of a pure ^D2p amplitude. This is precisely what 
one would expect in a classical picture. On the other hand in the central 
bump the a-particle distribution is isotropic and it is the pion distribution 
with respect to the beam direction which inherits this shape. 

The final expression for the differential cross section in the laboratory is 



dV 

dQ dpa 



1 plk* 



■ J dn*Y,\M*\^, (7) 



192(27r)5mdMxPd^a 

where for quantities denoted by an asterix (*) the pion momenta and angles 
are evaluated in the vrvr rest frame. Though we have only calculated tt+tt" 
production, from isospin arguments the corresponding cross section for neu- 
tral pions should be half of this. However, because of the narrowness of the 
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ABC peak, it is important to evaluate the kinematics separately in the two 
cases. 

In the numerical evaluation we use the S'-state Paris wave function |^ 
for the deuterons and the a:dd cluster wave function of Forest et al. iTD 



According to the latter, the I = 0, S = 1 S'-state nucleon-nucleon distri- 
butions look very similar to those of the deuteron when calculated with the 
corresponding nucleon-nucleon force provided that the separation r^p < 2 fm. 
The scale factor between the two intensities is 4.7. Since pion absorption in 
the deuteron occurs primarily when the nucleons are close together, we need 
only the number of such 'small' deuterons in the a-particle. Thus, assum- 
ing that the 'deuterons' in the a-particle have independent distributions, we 
normalise the Q;:dd wave function to A*"^ = 2.3. 

The only beam energy for which the Saclay group measured an angular 
distribution is Td = 1250 MeV and this corresponds to a region where 
the ^D2p transition pp — ^ d7r+ amplitude is dominant [Q. When compar- 
ing with such a-particle momentum spectra, it is important to take into 
account the momentum resolution and this has been done by smearing the 
theoretical predictions over a resolution function with a = 10 MeV/c. After 
doing this, our predictions are typically a factor of two too high and, for 
ease of comparison, they have been divided by that in fig. 0. Such a small 
reduction in overall normalisation could easily be due to initial distortion in 
the deuteron-deuteron system combined with uncertainty in the a-particle 
normalisation. 

The overall agreement between theory and experiment in fig. ^ is remark- 
ably good and suggests strongly that, at least in this channel, the ABC is 
indeed a kinematic enhancement due to two independent p-wave pion pro- 
ductions. At a laboratory angle of 11°, the cm. angle in the ABC peak is 
about 45° and so the predictions are following the change of about a factor of 
2.5 in cm. differential cross section. The small deviations at the valleys and 
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the central bump at the larger angles may be due to having limited ourselves 
to just one amplitude. Away from 1250 MeV data the ^D2P dominance as- 
sumption is less good and, though we reproduce the other Saclay data to 
better than a factor of two, the details of the momentum spectra must wait 
until we have included a more complete set of pp — ^ dvr"'' amplitudes in the 
calculation. This is currently in progress. 

A fuller set of input amplitudes is also required in order to predict the 
deuteron analysing powers which have been measured near the forward di- 
rection in our energy domain at the SPESIII spectrometer at Saclay . In 
particular the ei ■ 62 form in eq. (|l|) predicts zero deuteron tensor analysing 
power whereas experiment yields a small value but with significant structure 
which might arise through interferences with other pp d7r+ amplitudes. 

Extra confirmation of our approach might be provided by an exclusive 
measurement of the dd an^Ti^ reaction, since in the central region we 
predict that the angular distribution should behave like (3cos^6',r + 1)^- 
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Figure 2: The a-particle momentum spectra from the dd aX reaction for 
six different laboratory angles at = 1250 MeV. The experimental data 
of ref. H are compared with the predictions of our double-pion-production 
model after division by an overall factor of two. The broken lines are the raw 
calculations and yield the solid lines after smearing over the experimental 
momentum resolution. 
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